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From  this  analysis,  one  has  to  solve  a  system  of  ODE’S  for  tne 
parameters  of  the  soliton. 


This  has  been  done  on  the  VAX.  Vie  have  also  integrated  numerically  the 
full  problem  to  see  the  effects  of  radiation. 
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(T  +  Y .  ) W  =  vw  -  (fir  +  Si)W'W 
t  r  1  xx 

The  boundary  conditions  are  periodic. 

Different  parameter  values  are  used  to  check  the  stability  of  the  x- 
inaependent  and  the  soliton  solutions.  The  equation  is  integrated  in  time 
using  a  pseudo-spectral  method.  The  linear  part  of  the  equation  is  solved 
for  half  time  step  using  the  FFT  (the  Fast  Fourier  Transform).  The 
nonlinear  part  is  solved  exactly  for  the  remaining  half-time  step.  Data  is 
collected  and  analyzed. 
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Current  Research 

There  are  two  major  projects  that  form  a  focus  for  my  research  at 
the  moment.  The  first  is  a  stability  problem  for  travelling  waves  in  a 
fast-slow  system  of  reaction-diffusion  equations.  The  travelling  wave  is 
constructed  by  piecing  together  solutions  of  reduced  problems  and  the 
question  is  to  understand  how  the  stability  information  for  these  reduced 
problems  is  put  together  to  determine  the  stability  of  the  full  problem. 

This  is  related  to  my  earlier  work  on  the  FitzHugh-Nagumo  equations  but 
significantly  more  difficult  as  the  underlying  phase  space  dimension  is 
4  rather  than  3.  R.  Gardner  of  the  University  of  Massachusetts  proved 
the  existence  of  the  waves  and  the  extra  space  dimension  translates  into 
non-trivial  behavior  in  the  slow  part  of  the  system.  The  understanding 
of  the  effect  of  this  new  feature  on  the  stability  problem  is  quite 
subtle.  Gardner  and  I  are  working  on  this  problem  jointly. 

The  second  problem  relates  to  standing  wave  solutions  of  the  nonlinear 
Schrodinger  equations.  I  have  proved  an  instability  result  and  applied 
it  to  nonlinear  optical  waveguides  with  Moloney.  I  am  currently  involved 
in  extending  the  results  to  a  wider  class  of  waves.  But  what  is  more 
interesting  is  that  the  results  have  suggested  some  striking  connections 
between  the  methods  of  proof  used  to  find  the  waves  in  the  first  place, 
namely,  the  dynamical  systems  approach  (as  developed  by  Kiipper  and  myself) 
and  the  variational  approach  (as  develcped  by  Strauss  and  Berestycki-Lions ) . 
These  techniques  have  up  until  now  developed  independently  and  a  relationship 
between  them  will,  I  believe,  open  up  the  possibility  of  exciting  new  results. 


Instability  of  Standing  Waves  for 
Nonlinear  Schrodinger  Type  Equations 


by 


Christopher  K.R.T.  Jones 

Division  of  Applied  Mathematics 
Brown  University 
Providence,  Rhode  Island  02912 


Current  Address: 

Department  of  Mathematics 
University  of  Maryland 
College  Park,  Maryland  20742 


This  work  is  dedicated  to  the  memory  of  Charles  Conley 


This  research  was  supported  in  part  by  the  National  Science  Foundation  NSF 
Grants  DMS  8501961  and  DMS  8507056,  Air  Force  Office  of  Scientific 
Research  AFOSR  Grant  #83-0027,  and  Army  Research  Office  ARO  Contract 
DAAG-29-83-K-0029. 


Instability  of  Standing  Waves  in  Nonlinear  Optical  Waveguides 


C.K.R.T.  Jones 


Department  of  Mathematics 
University  of  Arizona 
Tucson,  Arizona  85721 


U.S.A. 


J.V.  Moloney 
Physics  Department 
Heriot-Watt  University 
Riccarton,  Edinburgh  EH14  4AS 


Scotland,  U.K. 


PACS  numbers:  02.78  42.20 


Abstract 

A  new  mathematical  instability  technique  is  presented  and  applied  to 
determine  the  stability  properties  of  a  physically  important  class  of 
standing  waves  in  nonlinear  planar  optical  waveguides.  The  method  is 
illustrated  by  a  case  where  soliton  perturbation  techniques  or  variational 
methods  are  inapplicable. 
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New  Class  of  Instabilities  in  Passive  Optical  Cavities 

D.  \V.  McLaughlin,  J.  V.  Moloney,  and  A.  C.  Newell 
Applied  Mathematics  Program  and  Optical  Science  Center,  University  of  Arizona,  Tucson,  Arizona  85  721 

(Received  2  July  1984) 

In  this  Letter  we  show  that  the  fixed  points  of  the  Ikeda  map  are  more  unstable  to  perturbations 
with  a  short-scale  transverse  structure  than  to  plane-wave  perturbations.  We  correctly  predict  ihe 
most  unstable  wavelength,  the  critical  intensity,  and  the  growth  rates  ol  tnese  disturbances.  Our 
result  establishes  that,  for  a  large  class  of  nonlinear  waves,  spatial  structure  is  inevitable  and  drasti¬ 
cally  alters  the  route  to  chaos.  In  an  optical  cavity  the  consequence  is  that  the  period-doubling  cas¬ 
cade  is  an  unlikely  scenario  for  transition  to  optical  chaos. 

PACS  numbers:  42.65.-k 

In  this  Letter,  we  announce  a  new  and  unexpected  in  which  diffraction  effects  are  neglected.  The  purpose 

result,  an  instability  whose  consequences  have  ramifi-  of  this  Letter  is  to  point  out  that  this  assumption  is  not 


cations  for  a  large  class  of  nonlinear  wave  problems 
whose  dynamics  can  be  described  by  envelope  equa¬ 
tions.  Specifically,  it  deals  with  periodically  forced 
field  equations  of  the  universal  nonlinear  Schrodinger 
type.  This  instability  changes  the  whole  character  of 
the  route  of  the  system  from  a  simple  to  a  turbulent 
state.  It  generates  spatial  structure,  and  the  subse¬ 
quent  onset  of  chaotic  behavior  completely  bypasses 
the  period-doubling  scenario  which  is  relevant  if  spa¬ 
tial  structure  is  ignored.  Moreover,  the  scenario  which 
does  emerge  has  a  universal  character  of  its  own.  Ex¬ 
amples  of  this  phenomena  are  found  in  optics,  either 
in  the  transmission  along  optical  fibers  or  in  optically 
bistable  cavities.1  It  is  in  the  latter  context  that  this 
work  is  presented. 

In  this  problem  we  are  interested  in  the  long-time 
state  of  a  continuous  laser  signal  which  is  recirculated 
through  a  nonlinear  medium  In  examining  one  par¬ 
ticular  manifestation  of  optical  bistability  (a  ring  cavity 
with  Kerr  nonlinearity),  Ikeda2  wrote  a  map  expressing 
the  (complex)  amplitude  g„  +  i  of  the  electric  field  E 
on  the  ( n  +  1  )st  pass  through  the  cavity  as  a  function 
of  electric  field  amplitude  on  the  nth  pass; 

£„  +  i  =  a  4-  «gnexp[/<&0+  ipLNU)/ 21.  (1) 

In  (1),  a  is  the  amplitude  of  the  input  field,  R  <  1  the 
reflectivity  of  the  mirrors,  <t> o  the  detuning  parameter, 
p  is  (effectively)  the  length  of  the  nonlinear  medium, 
and  /V  (&,£„*)  measures  its  nonlinear  response.  Two 
cases  are  usually  studied:  (1)  the  saturable  medium, 
Nil)  —  (1  +  2/)-1,  (2)  the  Kerr  medium.  Nil) 
=  —  1-4-2/,  which  is  the  small  intensity  limit  of  the 
saturable  case.  Equation  (1),  called  the  Ikeda  map,  is 
a  two-dimensional  invertible  map  and  exhibits  a 
variety  of  behavior  which  is  already  well  documented 
in  the  literature.2"4  In  various  parameter  ranges  (the 
two  parameters  which  are  varied  are  a  and  p),  one 
finds  multiple  fixed  points  (see  Fig.  1)  and  sequences 
of  period-doubling  bifurcations  leading  to  chaotic  at¬ 
tractors. 

The  map  (1)  invokes  the  plane-wave  approximation 


justified  even  for  cases  in  which  the  input  beam  is  very 
slowly  varying  in  the  transverse  direction  xix.y)  and 
the  Fresnel  number  is  large.  The  reason  is  that  the 
fixed  point  solutions  of  the  plane-wave  map  (1)  are 
more  unstable  to  perturbations  with  a  short-scale 
transverse  structure  than  they  are  to  perturbations 
with  plane-wave  structure.  To  emphasize  this  point, 
the  numerical  experiment  discussed  in  this  Letter  is 
run  at  a  parameter  value  p  for  which  the  fixed  points 
of  the  Ikeda  map  are  stable! 

This  discovery  has  important  ramifications.  First,  it 
shows  that  the  initial  bifurcation  of  the  system  intro¬ 
duces  an  extra  dimension  into  the  problem,  a  short¬ 
wave  transverse  excitation  of  temporal  period  two. 
This  extra  dimension  affects  significantly  the  subse¬ 
quent  behavior  of  the  system.  As  the  stress  parame¬ 
ters  are  raised,  no  period-doubling  cascades  into  chaos 


FIG.  1.  (a)  The  multivalued  response  of  the  amplitude  of 
the  Rxed  point  |g|  vs  a  at  fixed  values  of  p  for  the  Ikeda 
map.  (b)  A  graph  of  h  ip,  t  )  vs\/r  =  -JyK  for  p  —  pi  =  pgg’ 
equal  to  0. 1 1  and  0.24. 
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The  stability  of  stationary,  TE0-type,  nonlinear,  thin-film  guided  waves  was  investigated 
numerically  for  both  symmetric  and  asymmetric  planar  waveguides  with  nonlinear  cladding  and 
substrate  layers.  It  is  found  that  large  regions  of  the  dispersion  curves  are  unstable  at  high  powers. 


Unique  properties1-7  have  been  predicted  for  waves 
guided  by  thin  films  when  one  or  more  of  the  guiding  media 
exhibit  a  field-dependent  refractive  index.  Self-focusing 
bounding  media  lead  to  multiple  new  branches  with  power 
thresholds,  as  well  as  field  distributions  whose  maxima  shift 
from  the  film  to  the  bounding  media  with  increasing  pow¬ 
er. 1-7  In  fact,  this  geometry  has  been  identified  as  an  excel¬ 
lent  candidate  for  all-optical  switching,  with  our  without 
bistability.5-7  To  date,  however,  theoretical  analysis  has  been 
based  solely  on  steady-state  solutions  to  a  nonlinear  wave 
equation  which  contains  an  intensity-dependent  refractive 
index.  The  salient  question  is  whether  these  wave  solutions 
are  stable  on  propagation,  and  the  consequences  of  possible 
unstable  regions  to  proposed  devices.  For  the  related  prob¬ 
lem  of  self-focusing  of  plane  waves  in  infinite  media,  Kolo- 
kolov*  has  shown  that  the  solutions  are  stable  fordP/dP>0 
w  here  Pis  the  power  and/?  is  the  power-dependent  refractive 
index.  Numerical  propagation  studies  of  nonlinear  waves 
guided  by  the  interface  between  a  self-focusing  and  a  power- 
independent  medium  by  Akhmediev  and  co-workers9- l() 
have  led  to  a  similar  conclusion.  Recently,  a  theory  based  on 
phase  portraits  has  been  developed"  for  the  stability  of  TE„ 
waves  guided  by  a  thin  film  bounded  by  self-focusmg  media, 
and  the  important  conclusion  of  that  work  is  that  the  waves 
are  unstable  on  negatively  sloped  branches  [dP /d(3 <0)  of 
the  nonlinear  dispersion  curve,  and  that  they  are  stable  on 
positively  sloped  regions  provided  that  self-focusing  occurs 
in  only  one  nonlinear  medium.  In  this  letter  we  report  a  test 
of  this  conclusion  via  a  numerical  investigation  of  the  stabil¬ 
ity  of  TE„  solutions  for  films  bounded  by  self-focusing  me¬ 
dia. 

The  geometry  analyzed  consists  of  a  film  (\xl  <d,  re¬ 
fractive  index  n„)  bounded  by  two  nonlinear  media  with  low- 
power  indices  n ,  and  n:,  as  shown  in  Fig  1 .  Since  the  numeri¬ 
cal  analysis  is  performed  in  the  slowly  varying  phase  and 
amplitude  approximation,  we  write  the  optical  field  as 


E(r)  =  (IF(x^)e'1 


where /?  is  the  effective  guided  wave  index.  k0  =  to/c  and  the 
variation  in  the  amplitude  term  W(-x.r)  along  the  propaga¬ 
tion  direction  z  is  assumed  to  be  small  over  one  wavelength. 
The  refractive  index  in  the  various  media  is  given  by 


VS  .  »tt  48  (13)  3’ 


n:[  ]  W(xj)  ]-]  =  rr  -  a,. ;  W(xj)  : 

y  =1-2  ( jxj></) 


n2 [  j  ^'(x^)!2]  =  n-  (  x\<d)  .  (2b) 

Substituting  into  the  nonlinear  wave  equation  and  retaining 
only  the  first  derivative  of  H’(x^z)  with  respect  tor  leads  to 


dz  dx~ 


X{/?2  -n2[|H..xj:),:]}W'U,z)  =0. 


Analytical  stability  analysis  of  Eq.  (3)  is  complicated  by 
the  fact  that  it  is  a  partial  differential  equation  and,  further¬ 
more,  is  a  Hamiltonian  system.  The  usual  stability  analysis 
for  dissipative  dynamical  systems  does  not  apply  ( unless  one 
deliberately  introduces  losses  into  the  problem ) .  One  reason 
for  the  difficulty  in  studying  the  stability  of  Eq.  (3)  is  that 
many  of  the  eigenvalues  of  the  linearization  of  ( 3 )  lie  on  the 
imaginary  axis  which  is  precisely  the  condition  for  instability 
in  a  dissipative  dynamical  sy  stem.  Although  it  has  been  pos¬ 
sible"  to  perform  a  stability  analysis  for  TEm  with  m  —  0.  it 
is  necessary  to  proceed  numerically  for  m=*=  0.  First,  the 
steady-state  solutions  w  ith  W(xj)— -  lV(x )  were  obtained  in 
the  usual  way 1-7  to  ootain  a  field  solution  corresponding  to  a 
particular  point  on  one  of  the  nonlinear  guided  wave  solu¬ 
tion  branches.  This  distribution  was  then  assumed  to  be 
launched  at  z  =  0  and  Eq.  (3)  was  solved  numerically  for 


n?*a|E|? 


FIG  l  Nonlinear  *a\ej  .  iccrr.ctry  studied  for  stabiht> 
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Chaos  and  Turbulence;  is  there  a 
connection? 


by 
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Tucson,  Arizona  85721 


Abstract .  In  this  essay  we  discuss  the  relation  of  chaos,  which  is 
the  unpredictable  behavior  associated  with  finite  systems  of  ordinary 
differential  and  difference  equations,  and  turbulence,  which  is  the 
unpredictable  behavior  of  solutions  of  infinite  dimensional,  nonlinear 
partial  differential  equations.  The  evidence  that  there  is  some 
connection,  at  least  in  certain  regimes  of  parameter  space,  is 
sufficiently  convincing  to  provide  the  motivation  to  search  for 
analytical  means  for  reducing  the  governing  partial  differential 
equations  to  either  a  finite  system  of  ordinary  differential  equations 
or  a  much  simpler  partial  differential  equation  of  universal  type.  A 
successful  reduction  scheme  must  capture  the  spatial  structure  of  the 
dominant  modes  accurately  and  we  suggest  ways  of  finding  these 
structures  in  certain  limiting  situations.  Five  such  schemes  are 
presented  and,  in  each  case,  the  approximation  is  related  in  some  way 
to  the  presence  of  a  small  parameter,  near  critical,  nearly 
integrable  or  nearly  periodic.  One  of  these  reductions  leads  to  the 
complex  Ginzburg-Landau  equation,  which  has  universal  character,  and 
its  importance  is  stressed.  In  connection  with  this  equation,  we 
introduce  the  terms  "wimpy"  turbulence  and  "macho"  turbulence  to 
connote  the  crucial  differences  between  the  behavior  of  its  solution 
in  one  and  two  space  dimensions,  a  difference  which  has  much  in  common 
with  the  contrast  between  two  and  three  dimensional  high  Reynolds 
number  flows  because  of  vorticity  production.  In  the  final  section, 
several  ideas  concerning  the  nature  of  high  Reynolds  number,  fully 
developed  turbulence  are  presented  and  the  possible  roles  of  singular 
solutions  and  "fuzzy"  attractors  are  discussed.  Throughout  the  essay, 
we  argue  that  before  much  new  progress  is  made,  one  has  to  understand 
the  onset  of  spatial  chaos,  that  is,  the  transition  from  a  spatially 
regular  state  (possibly  with  a  chaotic  temporal  behavior)  to  one  in 
which  the  spatial  power  spectrum  is  broadband.  This  question  is  a 
major  focus  of  our  present  research  program. 

This  contribution  is  dedicated  to  the  memory  of  Dick  DiPrima,  a 
good  friend  and  long  time  colleague  who  left  us  too  soon. 
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The  geometry  of  the  Hill  equation  and  of  the  Neumann  system 

By  N.  M.  K.RCOLANI  j  AND  H.  Flaschka 
Department  of  Mathematics,  The  University  oj  Arizona,  Tucson ,  Arizona  85721,  U.S.A. 


Let  there  be  given  a  finitergap  operator  L  =  d2/dx2  +  q  and  its  Baker  function  f(x,  p), 
which  is  analytic  for  p  on  a  certain  hyperelliptic  curve  C.  It  is  shown  that  a  sequence 
of  Backlund  transformations  maps  C  to  a  projective  space.  This  embedding  can  be 
interpreted  as  a  matrix  representation  of  the  Hill  equation  by  the  Neumann  system 
of  constrained  harmonic  oscillators.  The  image  curve,  C ,  lies  on  a  rational  ruled 
surface;  the  structure  of  this  surface  is  explained  by  use  of  ideas  due  to  Burchnall  & 
Chaundv  {Proc.  R.  Soc.  Lond.  A  118,  557-583  (1928)).  Baker  functions  and  Backlund 
transformations  are  then  used  to  define  a  (many-to-many)  correspondence  between 
effective  divisors  on  the  curve  C  and  points  lying  on  a  quadric,  or  in  the  intersection 
of  two  or  more  quadrics.  This  relates  the  theory  of  the  Hill  equation  to  earlier  work 
of  Knorrer,  Moser  and  Reid.  It  is  then  shown  that  the  Kummer  image  of  the  Jacobian 
of  C  can  be  realized  as  a  hypersurface  in  the  space  of  momentum  variables  of  the 
Neumann  system.  Further  projects,  such  as  extensions  to  non-hypcrelliptic  curves,  are 
outlined. 


1.  Introduction 

A  differential  operator  L  =  D2 +  <?(*),  D  =  d/dx,  is  said  to  be  ‘finite-gap’  if  it  commutes  with 
a  differential  operator  B  of  odd  order,  (L,  B]  =  0.  A  ‘finite-gap  potential’  q(x)  is  therefore  a 
time-independent,  or  stationary,  solution  of  an  equation  8L/dl  =  [B,  L]  in  the  Korteweg-de 
Vries  hierarchy.  Because  L  and  B  commute,  they  have  a  common  eigenfunction 

L  ij/  =  Eif, 

Bif  =  Rif. 

The  eigenvalues  E,  R  are  known  to  be  related  by  an  algebraic  equation 

Ri  =  Vl‘  (E-cf,  (1) 

and  the  common  eigenfunction  if  (the  ‘ Baker  function’)  is  an  analytic  function  on  the  Riemann 
surface  (1)  or,  equivalently,  a  holomorphic  section  of  a  certain  line  bundle  on  (1).  Until  now, 
the  theory  of  finite-gap  operators  has  drawn  mostly  on  the  analytical  aspects  of  Riemann 
surfaces  and  on  their  abstract,  intrinsic  geometry. 

Our  aim  in  this  paper  is  to  explain  some  of  the  extrinsic  properties  of  the  curves,  line  bundles, 
and  isospectral  tori  (Jacobians)  when  those  are  embedded  as  concrete  objects  in  a  projective 


There  are  several  reasons  for  studying  geometric  realizations  of  the  finite-gap  operator  theory. 
The  classical  theory  ofeurves  and  Jacobians  is  very  beautiful,  and  an  interesting  statement  about 
abstract  curves  and  line  bundles  should  be  worth  repealing  about  concrete  representations. 
Furthermore,  when  one  integrable  system,  like  the  stationary  Lax  equation  [L,  B)  =0,  is 
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li  is  shown  that  the  structure  of  the  stationarv  dislocations  which  occur  in  a  wide  varietv  of  patterns  in  rmnequilibnunt 
s\ stems  is  given  b\  a  self-similar  solution  of  the  universal  Cross- Newell  equation. 


There  has  been  a  tremendous  resurgence  of  interest 
lately  in  the  structure  and  properties  of  symmetry 
breaking  instabilities  in  driven  systems  far  from  equi¬ 
librium.  The  planforms  and  patterns  which  emerge 
from  these  instabilities  are  observed  in  a  wide  variety 
of  physical  situations,  from  convection  in  fluids,  in 
liquid  crystals  and  in  binary  mixtures  to  the  solidifi¬ 
cation  processes  in  undercooled  liquids.  Of  particular 
interest  are  patterns  which  form  in  large  aspect  ratio 
systems;  for  example,  the  Rayleigh  Btfnard  convec¬ 
tion  of  a  Boussinesq  fluid  in  a  box  wide  enough  ro 
contain  many  rolls.  Although  for  certain  parameter 
ranges  and  in  infinite  horizontal  geometries,  there 
exist  stable,  fully  nonlinear,  spatially  periodic  solu¬ 
tions  of  the  Oberbeck- Boussinesq  equations 

u(x,t)=f(0,  A.R),  (1) 

which  correspond  to  a  field  of  infinitely  long,  straight, 
parallel  rolls,  the  patterns  which  are  typically  seen 
are  much  more  complicated,  involving  patches  of 
curved  rolls,  defects  such  as  roll  dislocations  and 
grain  boundaries.  In  an  effort  to  treat  the  statics  and 
slow  dynamics  of  these  patterns.  Cross  and  Newell 
(1 1  developed  a  theory  which  averages  over  the  de¬ 
tailed  local  structure  of  the  rolls  and  concentrates  on 
the  global  and  universal  properties  of  the  pattern  it¬ 
self.  The  starting  point  of  the  theory  is  the  existence 
and  stability  of  the  solution  ( 1 )  representing  the 
underlying  roll  structure.  Here/ is  2n  periodic  in  the 
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phase  0  whose  gradient  is  the  constant  wavevcctor 
k.  The  periodicity  demand  gives  rise  to  a  relation 
(equivalent  to  the  frequency  dependency  on  ampli¬ 
tude  in  nonlinear  oscillators) 

o(k,A.R)=  0  „  (2) 

between  the  wavenumber  k  =  |Ac|.  the  amplitude  A 
and  the  stress  parameter  R  which  in  the  case  of  con¬ 
vection  is  the  Rayleigh  number.  The  appearance  of 
the  modulus  of  k  reflects  the  rotational  symmetry. 

The  regions  of  stability  of  iJjc.se  solutions  in  the  R.  k 
plane  have  been  mapped  out  by  Busse  [2|  (the  Busse 
balloon). 

The  Cross  Newell  theory  develops  a  universal 
phase  equation  for  the  slowly  changing  wavevcctor  k , 
a  change  necessitated  by  the  influence  of  distant,  but 
finitely  distant,  horizontal  boundaries  to  which  the 
roll  axis  is  perpendicular  (if  the  thermal  contact  be¬ 
tween  wall  and  fluid  is  good).  For  order  one  values 
of  R  -  A'c  (Rc  is  the  transition  value  at  which  rolls 
first  appear),  the  amplitude  is  still  determined  algebrai¬ 
cally  in  terms  of  the  phase  gradient  by  (2).  In  the  ab¬ 
sence  of  the  mean  drift  effect,  which  we  do  not  con¬ 
sider  here,  the  phase  equation  is 

r(k)0t  +V(kB(k))  +  F(k)(D]D2  +£>2Dl >A 

+  ...  =  0.  (3) 

where  the  functions  r(k)  >  0,  R(k)  and  F(k )  arc  cal¬ 
culated  easily  from  the  particular  model  of  interest 
(see  ref.  1 1 1 )  and  the  operators  D  |  and  /T,  ate  2k- V 
+  V-Ar  and  V2,  respectively.  The  function  B{k)  always 
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CHAPTER  4 

AN  APPLICATION:  THE  ROLL-UP  OF  A  FINITE  VORTEX  LAYER 

The  evolution  and  regularity  of  vorticity  distributions  with 
varying  degrees  of  smoothness  are  of  great  interest  to  those  who 
study  the  Euler  equations  of  fluid  flow.  A  periodically  perturbed 
vortex  sheet,  for  example,  is  now  widely  believed  to  acquire  a 
curvature  singularity  in  finite  time.  Using  the  methods  developed  in 
Chapters  2  and  3  for  the  highly  accurate  application  of  boundary 
inteyral  methods  to  multiply  connected  domains,  the  regularity  of  a 
thin,  periodic  layer  of  constant  vorticity  is  investigated 
numerically.  Numerical  results  suggest  that,  like  the  vortex  sheet, 
the  interfaces  develop  a  curvature  singularity,  but  now  only  in 
infinite  time. 

I.  Background 

The  simplest  model  of  a  high  Reynolds  number  shear  layer  is  a 
surface  of  discontinuity,  or  vortex  sheet,  between  two  shearing, 
inviscid,  irrotational  fluids  as  pictured  schematically  in  Figure  11. 
The  instability  of  such  an  interface  is  prototypic  and  is  the  well- 
known  Kel vin-Helmholtz  instability.  For  simplicity,  assume  the  fluid 
is  of  uniform  density,  and  that  u  +Uj  as  y  -*•  +°°.  Then  the 

linear  evolution  of  a  normal  mode  of  amplitude  A(k,t)  and  wave 
number  k  is  given  by 
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This  paper  addresses  the  possible  connections  between  chaos,  the  unpredictable  behavior  of  solutions  of  finite  dimensional 
s\ stems  of  ordinary  differential  and  difference  equations  and  turbulence,  the  unpredictable  behavior  of  solutions  of  partial 
differential  equations.  It  is  dedicated  to  Martin  Kruskal  on  the  occasion  of  his  60th  birthday 


1.  Introduction 

The  chaos  that  occurs  in  p.d.e.  ’s 
cannot  be  fathomed  by  legalese 
so  n  r  apply  Occam ’s  razor 
and  using  a  laser 
study  structures  in  ring  cacities 
An  appealing  idea  of  modern  dynamics  is  that 
the  complicated  and  apparently  stochastic  time 
behavior  of  large  and  even  infinite-dimensional 
nonlinear  systems  is  in  fact  a  manifestation  of  a 
deterministic  flow  on  a  low-dimensional  chaotic 
attractor.  If  the  system  is  indeed  low  dimensional, 
it  is  natural  to  ask  whether  one  can  identify  the 
physical  characteristics  such  as  the  spatial  struc¬ 
ture  of  those  few  active  modes  which  dominate  the 
dynamics.  Our  thesis  is  that  these  modes  are  closely 
related  to  and  best  described  in  terms  of  a.  vmptot- 
ically  robust,  multiparameter  solutions  of  the  non¬ 
linear  governing  equations.  We  find  it  hard  to 
define  this  robust  nature  precisely,  but  loosely 
speaking  the  idea  is  that  these  solutions  are  very 
stable  and  resilient  asymptotic  states.  They  may  be 
coherent  lumps  like  solitons  and  solitary  waves. 
They  may  have  the  form  of  coherent  wave  packets. 
They  may  have  self-similar  form.  They  need  not 
necessarily  be  the  asymptotic  states  which  develop 
as  /  tends  to  infinity;  structures  which  develop 
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singularities  in  finite  time  like  those  involved  in  the 
collapse  of  Langmuir  waves  or  in  filamentation  in 
nonlinear  optics  are  also  candidates.  For  example, 
singular  solutions  of  the  Euler  equations  may  be 
useful  in  understanding  the  behavior  of  the 
Navier-Stokes  equations  at  high  Reynolds  num¬ 
bers.  Singular  solutions  like  defects  and  disloca¬ 
tions  certainly  do  play  important  roles  in  the 
pattern  formations  arising  in  continuum  and  con¬ 
densed  matter  physics.  The  key  idea  is  that  each  of 
these  structures  is  a  natural  asymptotic  state  that, 
by  virtue  of  the  various  force  balances  in  the 
governing  equations,  develops  an  identity  which 
docs  not  easily  decay  or  disperse  away. 

One  can  envision  two  types  of  chaos  occurring. 
The  first  is  a  phase  or  weak  turbulence  which 
arises  when  there  is  an  endless  competition  be¬ 
tween  equally  resilient,  localized  coherent  struc¬ 
tures  which  are  infinite  time  asymptotic  states  and 
which  are  initiated  at  random  at  various  parts  of 
the  physical  domain.  Examples  of  this  type  of 
turbulence  are  solitary  waves  in  the  one-dimen¬ 
sional  complex  envelope  equation,  Rayleigh- 
Benard  roll  patterns  with  different  orientations 
and  the  oscillatory  skew  varicose  states  in  low 
Prandtl  number  convection.  It  is  to  be  expected 
that  such  dynamics  may  be  low  dimensional.  The 
second  type  of  chaos  is  much  more  dramatic  and. 
for  want  of  a  better  word,  may  be  described  as  an 


0167-27H9/86/S03.50  *  Elsevier  Science  Publishers  B.V. 
(North-Holland  Physics  Publishing  Division) 


